Power law behavior for the zigzag transition in a Yukawa cluster 



We provide direct experimental evidence that the one-dimensional (ID) to two-dimensional (2D) 
zigzag transition in a Yukawa cluster exhibits power law behavior. Configurations of a six-particle 
dusty (complex) plasma confined in a biharmonic potential well are characterized as the well 
anisotropy is reduced. When the anisotropy is large the particles are in a ID straight line config- 
uration. As the anisotropy is decreased the cluster undergoes a zigzag transition to a 2D configu- 
ration. The measured dependence of cluster width on anisotropy is well described by a power law. 
A second transition from the zigzag to an elliptical configuration is also observed. The results are 
in very good agreement with a model for particles interacting through a Yukawa potential. 
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Plasma is a quasi-neutral gas that contains electrons, ions and (often) neutral atoms and 
which exhibits collective behavior [l| . A dusty (complex) plasma is a systems of interacting 
microscopic dust particles immersed in a standard electron-ion plasma. For typical labora- 
tory conditions, dust particles acquire a net charge g < from the electron and ion currents. 
The Coulomb interaction between the charged dust particles is shielded by the response of 
the free charge in the plasma. In the plane of a two-dimensional cluster, the inter-particle 
potential is well approximated by a Yukawa (i.e., Debye or shielded Coulomb) potential ^ 

V(r) = ^^e-''/^^, (1) 
^ ' 47reo r ' ^ ^ 

where A^i is the Debye screening length. An unshielded Coulomb interaction is recovered 
as Az) — 7- oo. Dusty plasma is an exceptional experimental system for examining the static 
and dynamic properties of Yukawa clusters-small clusters of particles interacting through a 
Yukawa potential. 

In laboratory experiments the dusty plasma is usually confined, typically through a 
combination of gravitational and electrical forces. Monodisperse microspheres may then 



or even three- 



form strongly-coupled one-dimensional (ID) [sHSj, two-dimensional (2D) [( 
dimensional systems j^. The general 2D potential well is biharmonic for small displacements 
from the potential energy minimum A dusty plasma confined in a 2D biharmonic 

well may be either in a ID straight line configuration, or a 2D configuration, where the 



transition from ID to 2D states is via the zigzag 



instabi 



12 



ity mm 



m]i5\. 



by assuming n identical particles 



A cold 2D dusty plasma can be modeled |10 . 
with mass m and charge q, where the ith particle is at {xi,yi). The potential energy [/ of a 
configuration {xi,yi} is given by the sum of the particles' potential energy with respect to 



the biharmonic well plus the sum of the potential energy for each pair-wise interaction, 



U 



i=l 



i=l j>i 



(2) 



lij 



where r 



;j — Xj)^ + {ni ~ UjY is the separation between particles i and j, and and 
Uy are the single-particle oscillation frequencies in the x and y directions, respectively. The 
equilibrium configuration minimizes U. 

We nondimensionalize Eq. using the variables = Xj/rg, rji = yi/vQ, pij = Vij/rQ and 
parameters 



Here the characteristic length 

is defined using the longitudinal frequency w^.. The dimensionless potential energy is then 

1=1 1=1 j>i ■' 

where Uq is the characteristic potential energy. Equilibrium configurations {rji,^i\ of the 
model [Eq. (E])] depend on three parameters: ra, /t and a^. Here k is the Debye shielding 
parameter and k = corresponds to an unshielded Coulomb interaction. The well anisotropy 
parameter is a^. We assume Uy > uj^ oi > 1 . Straight line configurations Vi = rji = then 
lie in the x (longitudinal) direction and are independent of a^. An equilibrium configuration 
{^i, rji} can be compared to an experimentally measured configuration {xj, i/i} by multiplying 
each coordinate by tq. 



Within the context of this model [Eq. (j5])], abrupt ID to 2D transitions 10|, uM 



predicted to occur for critical values of the three parameters n, k and due to a zigzag 



instability. This has been confirmed experimentally as n was varied by Melzer and 
Sheridan and Wells 16|. Melzer Uj also observed a zigzag transition as the neutral gas 
pressure was changed, which presumably changed the plasma density and therefore both k 



and a^. Sheridan and Wells 



16| confirmed that the zigzag transition caused by increasing 



n in small clusters behaves as a continuous phase transition between ID and 2D states, 
where the cluster width near the transition obeys a power law. Using the model of Eq. 
©. they also predicted Q that the .gzag transifonin s,„a.. clusters exh.bits power law 
behavior as either k or a is varied. Piacente, et al. [17] reached a similar conclusion from 
a computational study of unbounded Yukawa systems (n — ?■ oo) with a transverse parabolic 
well and longitudinal periodic boundaries. Structural transitions due to changes in k or 
have not been characterized experimentally. 

In this work we investigate the transitions a Yukawa cluster (i.e., a small dusty plasma) 
undergoes as the well anisotropy is varied. We use a novel experimental method in which 
the physical geometry of the confining potential well is changed without disrupting the 
plasma discharge. The discharge parameters are nearly constant, so that the Debye shielding 
parameter k is also nearly constant. The well anisotropy is accurately determined using 
measured center-of-mass frequencies for the dusty plasma. As predicted fl6^, the transverse 
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Figure 1: Schematic of the experimental geometry (not to scale). Six dust particles are confined in 
the biharmonic potential well created by a rectangular depression placed on the powered electrode 
in an rf discharge. The confining well anisotropy is varied by changing the width d while the plasma 
is on, and without losing the particles, using a lead screw. 

cluster width yrms exhibits power law behavior vs following the zigzag transition. The 
Debye shielding parameter is found by comparison to the model, and Debye shielding is 
shown to be important. 



The dusty plasma was studied in tlie^DO_?v 
imenT) apparatus (Fig. [1]) j?, Is. Il2 



16 
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UT (Dusty Ohio Northern University exper- 



19| . Argon gas was leaked into the vacuum 



chamber and the pressure was stabilized at 15.4 mtorr. At this pressure normal modes of the 
dust clusters are under damped [ist. Inside the chamber a capacitively-coupled rf electrode 
is used both to sustain the discharge and to levitate the dust particles. The rf frequency 
was 13.56 MHz, the forward rf power was ~ 7 W and the dc self-bias on the electrode was 
-90 V. 

The confining well (Fig. [1]) is produced by a rectangular aperture placed on the 89-mm 
diameter aluminum electrode, creating a biharmonic well in the horizontal (x-y) plane j^. 
The aperture is made of four aluminum bars, each with a cross section of 6.35 mm x 12.7 mm. 
The depth of the trap is 6.35 mm, the length is fixed at 50.8 mm, and the width d is adjustable 
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from to 38 mm. Here d can be adjusted while the plasma is on, and without losing the 
dust particles, by mean of a lead screw which is rotated from outside the vacuum. The lead 
screw is made of joined left-handed and right-handed threaded rods [unified national coarse 
(UNC) thread size 6-32] that move the two confining bars simultaneously closer together or 
farther apart a distance of 1/16 inch {k. 1.59 mm) per revolution, thus keeping the minimum 
of the potential well centered with respect to the disk electrode. This arrangement allows us 
to continuously vary the anisotropy parameter while maintaining nearly the same plasma 
discharge parameters. 

The dust particles used were melamine fornialdehyde resin spheres with a nominal diam- 
eter of 9.62 ± 0.09 /im. As noted previously [18], we believe that the particle diameter is 
actually 8.94 /im, giving m = 5.65 x 10"^^ kg. Confined dust particles are illuminated by a 
red diode laser sheet. A top- view camera mounted about 30 cm above the electrode is used 
to acquire dust particle position data. For this experiment, the resolution of the camera 
and associated optics was 16.47 /im/pixel. Sequences of 2048 frames were recorded at 27.9 
frames/s. A side-view camera was used to confirm that there were no out-of-plane particles. 
We acquired 20 data sets with the same n = 6 particles for trap widths d = 12.8 to 36.6 
mm. 

Representative measured configurations are shown in Fig. |2]for increasing d. As seen in 
Fig. El^a), the particles are in a straight line configuration when d is small, or equivalently, 

is large. As we increase d the cluster goes through two structural transitions, straight 
line to zigzag [Figs. Mjo) and (c)], and then zigzag to elliptical [Fig. EJ^d)]. The transition 
to the zigzag state occurs at ^ 23 mm and the transition to the elliptical state occurs at 
d ^ 30 mm. In the elliptical state all six particles lie on the convex hull of the cluster. 

To analyze the experimental data, we first find the dependence of the anisotropy param- 
eter on d. This is done by measurin g th e frequencies Ux and Uy of the cm. modes from 



the particles' Brownian motions jsl, Q, |16|. The power spectra of the time histories of the 
cm. coordinates, Xcm = (l/"-) Ucm = i^/n) were fitted with the expression 

for a driven damped harmonic oscillator to determine the corresponding normal mode fre- 
quencies Ux and Uy. For each d, the dust particle configuration was characterized by the 
root-mean-squared values Xrms and i/rms averaged over all frames, where 
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Figure 2: Measured particle configurations vs increasing trap width d, or equivalently, decreasing 
well anisotropy a^. (a) Straight line configuration, (b) and (c) zigzag configurations and (d) 
elliptical configuration. Values of the anisotropy parameter were determined from measured 
cm. frequencies. 

Here Xrms represents the longitudinal cluster size (i.e., its length), while i/rms is the trans- 
verse cluster width. In particular, yrms can be used as an order parameter for the zigzag 
transition sinceVrms = in the ID straight line configuration and y-rms > for the 2D zigzag 



configuration 



16| 



The measured cm. frequencies and Uy are plotted vs d in Fig. [3](a). Here (the 
frequency associated with the fixed 50.8 mm separation) shows a weak linear increase with 
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d, while ooy is roughly constant for (i < 20 mm and then decreases as d increases. That is, 
for these experimental conditions decreasing d below ~ 20 mm does not increase Uy because 
the sheath edge can no longer conform to the rectangular depression, but rather, is pushed 
upward This explanation is given further weight by measurements of the cluster height 
h above the top of the confining bars [Fig. El^a)]. Here h was determined by the height of the 
laser sheet used to illuminated the particles and has an uncertainty ±0.2 mm. As d increases 
h decreases, indicating that the sheath edge is moving into the rectangular depression. 

Rather than taking and Uy for each d and directly calculating a^, we fitted a line to 
and a quadratic curve to Uy for d > 20 mm, which spans the zigzag and elliptical transitions. 
The average value of over this range is 7.47 rad/s. By dividing the two fitted curves we 
determined as a function of d [Fig. E^b)]. We cover a wide range of anisotropics, where 

decreases with increasing d from ~ 14 to 2.5. If we could increase d to 50.8 mm to 
achieve an isotropic well, we would expect = 1. The slight increase of with d implies 
that K, oc Ux"^^^ increases weakly with a^. If the only change in k, is due to Ux, then we 
estimate that n increases by ^ 6% over the range of values considered. 

The measured cluster width i/rms and length Xrms are plotted vs in Figs. Hl^a) and (b), 
respectively. From the yrms data we observe a transition from the ID configuration to a 2D 
zigzag configuration at ~ 12. From ^ 12 to 6.5, the yrms data is well fitted by a power 
law with a critical value of. = 12.1 and an exponent of 0.35. The model solution shows the 
same power law behavior 16|]. The divergence from the power law fit at < 6.5 indicates 



a second transition 



16(1, from the zigzag configuration to an elliptical configuration. The 



elliptical structure 12| occurs when the second (fifth) particle in the cluster moves between 
the first and third (fourth and sixth) particles, thereby moving onto the convex hull of the 
cluster. Since the equilibrium shell structure for six particles in an isotropic well is (1, 5), a 
further symmetry breaking transition is needed to reach the o? = \ configuration. 

For the first value of d after the zigzag transition [o? = 11.8), the cluster was seen to flip 
between the two anti-symmetric zigzag conflgurations. (Since we can vary continuously 
it is possible to tune for this behavior.) A zigzag conflguration has two degenerate states 
since the energy is invariant under a reflection of particle coordinates y — )■ — In a stable 
zigzag conflguration these two states are separated by a potential barrier corresponding to 
the higher energy of the unstable straight line conflguration and the zigzag state is bistable. 
Just above the zigzag transition the barrier between the two states is low, and thermal 
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Figure 3: (a) Dependence of measured cm. frequencies uj^ and ujy and the height h of the cluster 
above the top of the weh on the width d. Here ujx is fitted with a hnear function and ujy with a 
quadratic curve, (b) Dependence of the smoothed anisotropy parameter = w^/w^ on d. 



fluctuations or small periodic perturbations can cause the system to flip between the two 
zigzag states, as was observed. For the next lower value of flipping was not seen, indicating 
that the potential barrier between the two degenerate zigzag states was then too high. 
Theoretical configurations which can be compared to the experiment are found by mini- 



mizing Eq. allowing us to predict Xrms = ^o^rms and i/r 



roVr 



function of a . 



To calculate a configuration three parameters, a^, n and n, are required. We know n, and 
has been measured, so we only need to determine the shielding parameter n. To do this 
we assume that k is independent of and then minimize the sum of the squared differences 
between the measured and computed values of yrms, giving k = 2.4 and vq = 1.25 mm. The 
Debye length A^) = tq/k, = 0.52 mm is comparable to the inter-particle distance and the 
average charge q = —1.2 x lO^e. These values are consistent with those found in previous 



experiments for similar plasma conditions 
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Figure 4: (a) Measured cluster width Urms (open circles), power law fit to data near the zigzag 
transition (dashed line) and widths computed from the model for k = 2.4 and rg = 1.25 mm (solid 
line) vs the measured anisotropy parameter a^. (b) Measured cluster length Xrms (open circles) 
and lengths computed from the model (solid line) vs a^. (Inset) Comparison of the percentage 
difference between the Xrms data and the model (open circles) with the percentage change in 
ro ~ Wx^^^ (broken line) vs a^. 

As shown in Fig. Hl^a), the cluster width Urms calculated using the model for k = 2.4 and 
tq = 1.25 mm exhibits excellent agreement with the experimental data, including reproduc- 
ing the transition between the zigzag and elliptical configurations. The agreement between 
the measured and predicted values of x^ms provides a cross check on the model. We find 
that the predicted values of Xrms [Fig. 11(b)] show good agreement for < 4, but system- 
atically underestimate the measured cluster length as increases. For the longest clusters, 
which are in a straight line configuration, the predicted cluster lengths are about 6% below 
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the measured values. As shown in the Fig. inset), this difference can be attributed to 
an experimental increase in tq ~ w^^^^ due to the measured decrease of Ux with a^. This 
implies that the average particle charge q is effectively constant, as previously observed for 



constant neutral pressure |8|, [l9|. This increase in tq with «^ is not as apparent for yrms 
since yrms ^ 0.2xrms and yrms — )■ as increases. 

In summary, we have provided direct evidence that the width of a Yukawa cluster exhibits 
power law behavior for the ID to 2D zigzag transition caused by decreasing the confining 
well anisotropy parameter a^, confirming a previous prediction [l^. Experiments were 
performed using a dusty plasma with n = 6 particles confined in the biharmonic well above 
a rectangular depression. The width d of the rectangular depression was increased while 
the plasma remained on to decrease while the Debye shielding parameter k remained 
essentially constant. The dependence of on d was accurately determined by measuring 
the cm. frequencies of the dusty plasma. A transition from the zigzag configuration to an 
elliptical configuration was also observed. The cluster width was found to be in excellent 
agreement with the predictions of a model which assumes identical particles confined in a 2D 
biharmonic well and interacting through a Yukawa potential. From the fit to the model we 
found the Debye length is comparable to the inter-particle distance, so that Debye shielding 
significantly effects the physics of these clusters. 
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